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EXPERIMENTAL ANALYSIS OF (MEP IN A ROTARY 
COMBUSTION ENGINE 

by H. J. Schock# W. J* Rico and P, R. Mong 

National Aoronouttcs and Space Administration 
Low^s Rosoareh Center 
Cleveland/ Ohio 44135 


SUMMARY 

This experimental work demonstrates the use 
of a NASA designed, real time Indicated Mean 
Effective Pressure (IMEP) measurement system 
which will be used to judge proposed improve- 
ments in cycle efficiency of a rotary combus- 
tion engine. This is the first self-contained 
instrument that is capable of making real time 
measurements of IMEP in a rotary engine. Pre- 
vious methods used require data recording and 
later processing using a digital computer. The 
unique features of this instrumentation include 
its ability to measure IMEP on a cycle by 
cycle, real time basis and the elimination of 
the need to differentiate the volume function 
in real time. 

Measurements at two engine speeds (2000 and 
3000 RPM) and a full range' of loads are pre- 
sented, although the instrument was designed to 
operate to speeds of 9000 RPM. 



THE ROTARY ENGINE is of interest as a 
future powerplant for general aviation air- 
craft because of its advantages in simplicity, 
compactness, high specific power, stratified 
charge adaptability and low noise levels. The 
main disadvantage of rotary engines is that 
they have a 12% higher specific fuel consump- 
tion (SFC) than do comnarable reciprocating 
engines according to work recently completed 
under a NASA contract (1)*. In an aircraft, 
the advantages of high specific power together 
with the multifuel tolerance of the rotary 
engine may outv/eigh the apparent disadvantage 
in SFC. 

To evaluate potential performance and 
economy gains that might be realized by fur- 
ther development of the rotary engine, a re- 
search program aimed at quantifying these 
possible gains was begun at NASA Lewis Re- 
search Center. It was decided uhat in order 
to evaluate potential gains and determine 
where losses occurred, basic measurement such 
as IMEP were necessary. 

The need to make IMEP measurements has 
long been recognized and much discussed. 

Brown (2,3) presented an extensive analysis on 
the measurements and interpretation of IMEP 
measurements in both direct injection and 
prechamber diesel engine Trown's approach 
uses two analog signals tiuau represent pres- 
sure and volume; analog circuits operate on 
these signals and produce an output voltage 
proportional to IMEP. Brown considered using 
digital circuitry but felt at the time that 
these circuits were not fast enough to give 
good detail and thus good accuracy. In later 
work, Lancaster et al (4) describes a proce- 
dure for measuring IMEP by averaging a digi- 
tized pressure signal at one crank angle. 

With this type of system, as with the method 
described by Brown, analysis of cycle-to-cycle 
variability and the causes of this variability 
are not possible as events within a cycle may 
be obscured by this averaging. 


*Numbers in parenthesis designate References 
at end of paper. 
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Researchers at NASA leRO have been making 
real time IMEP measurements in piston engines 
for several years. Rice and Cassidy (5) have 
been able to observe above mean values of IMEP 
in a cycle following a below mean value of 
IMEP. Tnis was attributed to the higher con- 
centration of unburned fuel in the residual 
gases from the first cycle, making the charge 
of the second cycle more fuel-rich. 

Attempts at identification of the causes 
of this cycle to cycle variability are under- 
way by many researchers. It is known that 
mixing and turbulence directly affect cycle to 
cycle variability, but the exact form o'* this 
relationship is complex and unknown. It will 
be especially important to quantify these 
variations as stratified charge engine devel- 
opment progresses. 

Several attempts have been made to measure 
IMEP in a rotary engine (^,7,8) but all have 
involved post run processing of collected 
data. One of the disadvantages to this method 
is that considerable periods of time may be 
spent collecting data that is defective before 
the researcher is aware of the erroneous re- 
sults. The objective of this investigation 
was to develop a system to measure IMEP on a 
cycle-by -cycle, real-time basis. The instru- 
ment could then be used to assess efficiency 
Improvements associated with changes in engine 
variables such as contustion chamber geometry, 
seal changes, and intake and exhaust manifold 
configuration without additional post-test 
data processing. In addition, because of its 
real time computational ability, this instru- 
ment could be used for engine control purposes. 

STATEMENT OF PROBLEM 

This experimental work was conducted to 
demonstrate the performance of an IMEP mea- 
surement system for use in a rotary coirbustion 
engine. Figure 1 is a graphic representation 
of how the rotor moves within the 
epitrochoidal shaped housing for a complete 
engine cycle. For the Wankel type rotary 
engine, a thermodynamic cycle consists of 1080 
crank angle degrees. These pictoral repre- 
sentations were produced via computer code 
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d«?veloped as part of this work. The equations 
which describe the housing profile and rotor 
motion were taken from Ansdale (9). As this 
fiq<.ii'e shnw«, in order In roUew fhr nrenMire 
in a thermodynamic system of interest througli- 
out the cycle, either a transducer must be 
mounted on the rotor face or o series of 
transducers located in the housing nxist be 
electronically switched to follow the rotor 
face throughout the cycle. The system chosen 
is the latter method and it uses four piezo- 
electric transducers located as shown. The 
unique feature of this instrumentation is its 
ability to compute and display IHEP on a cycle 
by cycle, real-time basis. 

EQUIRvlENT AND TEST ROCEDURE 

The test rig used for these experiments is 
shown in Figure 2. A 1978 Mazda 12A non- 
emission type engine rated at 82 KW (110 BHP) 

55 7000 RPNi was used. The two rotor engine has 
a displacement of I.IS liters (70 in^) and a 
9.A;1 compression ratio. This engine was 
chosen as it represents an inexpensive testbed 
for instrument developi'nent compared to avails 
able alternatives. In addition, it couples a 
late version of apex seal technology with an 
advanced contustion chamber design (10) to 
produce an engine with high specific power and 
good fuel economy relative to previous rotary 
engines. This engine is the same model as 
that used by Danieli (B).A four speed manual 
transmission was utilized to couple the engine 
to a 250 HP dynamometer. Motoring of the 
engine was accomplished with a 25" HP constant 
speed (1770 RPM), nx)tor located on the dynamo- 
meter. Hie engine was capable of being 
motored at several speeds by using various 
combinations of belt-driven pulleys between 
the motor and dynanxsmeter. Transmission 
gearing provided additional flexibility in 
choosing the desired engine speeds. 

Complete conventional instrumentation was 
installed on the engine in addition to the 
special combustion instrumentation which will 
be discussed in detail later in this report. 
Automatic control systems were utilized for 
engine speed, dynamometer load, and coolant 
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discharge temperature. Load was measured both 
with a reaction load cell and a Lebow Associ- 
ate Torquemeter model 160^-2K which was in- 
stalled between the dynamometer and the trans- 
mission as shown in Figure 2. An optical 
shaft angle encoder was installed on the front 
of the engine and was belt driven at a 3:2 
ratio by the eccentric shaft. 

Exhaust emissions were measured by a Scott 
Research Laboratories Test Console which in- 
corporated the following equipment: 

A Scott Total Hydrocarbon Analyzer Model 
A15 


A Scott Model 25Q 
A Scott Model 125 
A Beckman Model 864 
A Beckman Model 865 


O 2 Analyzer 
NO/NOx Analyzer 
CQ 2 Analyzer 
CO Analyzer 


Engine operation was monitored on the 
alpha-numeric CRT displays in the the control 
room. A micro -processor at the test cell was 
utilized to interface with a central minicom- 
puter which performed on line calculations and 
continuously updated the CRT displays. An 
IBM-3033 was used for necessary post-run pro- 
cessing and graphical analysis of collected 
data. 

The engine was operated at steady state 
conditions (speed, load and water tempera- 
ture), at which time the data were recorded 
and oscilloscope traces were photographed. 
Steady state was said to exist when there was 
agreement of _+ 5% between the air fuel (A/F) 
ratio from the direct measurement of air flow 
and fuel flow and the A/F ratio calculated 
from the method of Spindt. The dynamoireter 
'motor was then started and allowed to come up 
to the engine speed. At this point the fuel 
was cut off to the engine. After allowing a 
short time for stabilization, the nratoring 
data were recorded while the engine was still 
at normal operating temperature. 
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DEFINITION AND BASIC EQUATIONS 


In this work, the following definitions 
have been used. The subscript "F" indicates 
firing, indicates motoring and "0" indi- 
cates the quantity was back calculated from a 
torque anC speed measurement. 


IMEPT 



V/Vd 


s IMEP + PMEP 


IMEP - Indicated mean effective pressure 
(Fig. 3) 


PMEP - Pumping mean effective pressure 
(Fig. 3) 


FMEP - Friction mean effective pressure 


BMEP - Brake mean effective pressure 


MMEP ~ Motoring mean effective pressure 
BMP - Brake Horsepower 


The calculation of BMEP and MMEP for the 
rotary engine follows from the classical def- 
inition of brake mean effective pressure in 
any positive displacement engine. That is the 
theoretical constant pressure exerted during 
each powerstroke of the engine to produce the 
observed work output, assuming no losses. 

For a positive displacement rotary engine: 


Work 

* Strokes = Power Out 

(1) 

Stroke 

Min 


(B^EP 

* Vd) in-lbf * N/x = rhp 


(12 in/ft) (330Q0 ft-lbf) 

(2) 


min 


BMEP 

= (BHP) (12) (33000) X 

(3) 


Vd (n) 


Vd = Volume displaced during power 
stroke 

N s Speed of power shaft, RPM 
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X c Number of revolutions required for 
the work (Vq * BMfP) to be 
delivered . 

X 1 for this engine 

The calculation of MMEP follows directly. 

The following equality can now be written. 

BMEPd,F = 1MEP(t + PMEPp - FMEPp (4) 

PMEPf is a negative quantity for this natu- 
rally aspirated engine. 


Assuming FWSPp = FMEP^ 


FMEPp a MMEP 



( 5 ) 


The objective of motoring the engine is to 
determine FNEP^v^. The standardized work 
input, MMEP, consists of the following parts; 


PMEP|vi - PuiTping Work (neg term) 

FMEPM “ Engine Friction (neg term) 

QMSPM - Heat Transfer Losses (neg term) 

For the rotary engine EMEP is a significant 
factor in this calculation and is included in 
it as it represents work lost by means other 
than engine friction. 

One can then use the following equality to 
determine if this direct measurement corre- 
sponds to th's back calculated value. 



The instrument system described in this 
report can measure the following quantity in 
either the motoring or firing configuration. 


, PriV = IMEP + PMEP (7) 

Vd 


It can also be used to measure only the IMEP 
or only the PMEP in either the motoring or 
firing configuration, 
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REAL TIME IMEP MEASUREMENT SYSTEM 


PRESSURE TRANSDUCERS - This instrument 
system has been designed to be used with four 
transducers that are sequentially switched in 
synchronization with the rotor rotation. 

Since a rotor face spans a 120P sector, at 
least three transducers ore required, Because 
of geometric restrictions and the desire to 
have an overlap region, it is more practical 
to use four transducers. 

Piezoelectric transducers were used be- 
cause of their high frequency response and 
good high temperature performance. However, 
this type of transducer has no inherent refer- 
ence point and tends to have a high output 
offset drift with time and temperature. 

The ino'.vidual transducers are located in 
the rotor housing in such a manner that there 
is physical overlap by the rotor of two suc- 
cessive transducers (leading and trailing) for 
a brief period prior to switching. During 
this overlap period the average difference 
between the two transducers is measured and is 
used to offset the trailing transducer. Thus, 
when the correlator switches from the leading 
transducer to the trailing transducer a con- 
tinuous waveform is produced. The difference 
averaging circuitry is unique in that the 
averaging period is a function of the rotor 
angle' only and is thus independent of engine 
RPM. (iteration of the signal correlator ex- 
tends to 9000 engine RW (3000 rotor PPM), and 
a comprehensive description of its operation 
is given in Appendix A. 

The physical location of the pressure 
transducers is shown in Figure A. Transducers 
1, 2, and 3 are Kistler Model 601B and are 
cooled by the engine water jacket cooling 
system. Transducer number A is located in the 
exhaust manifold just downstream of the ex- 
haust exit on the rotor housing. Temperature 
of this transducer was maintained by a water- 
cooled copper wrap around the transducer 
holder. 

Flush mounted water cooled transducers 
would have been a more acceptable installa- 
tion, but the primary objective of this work 
was to demonstrate the feasibility of using 
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the system described, Geonnetrical constraints 
which would have made extensive rotor housing 
modification necessary did not seem appro- 
priate for proof of this concept. 

The sensitivity of piezoelectric trans- 
ducers changes with temperature, and thus it 
was necessary to determine their operating 
temperatures. To do this, four transducers 
were drilled #;nd fitted with thermocouples as 
shown in Figure 5. The thermocouples were 
placed as close as possible to the diaphragn 
of the transducer. Tests were then conducted 
at loads varying from 20 to 100 ft-lbf and 
speeds of 1000 to 5000 RPM to determine the 
temperature at which the transducers would 
operate , 

Based on the results of the temperature 
data each of the transducers installed had the 
sensitivity of its charge amplifier adjusted 
to compensate for the temperature at which it 
would operate. At 180° F coolant tempera- 
ture, the average temperature that the three 
transducers in the housing operated was 
193.10 F. The average temperature of the 
transducer in the exhaust manifold was 73.5° 

F. Spot checks of transducer temperatures 
while motoring the engine, which was done in 
each case immediately after firing, showed 
little deviance from the firing temperatures. 
All iDeasured temperatures are within the oper- 
ating limits prescribed by the transducer 
manufacturer (11). 

SIGNAL CORRELATOR AND SHAFT ANGLE 
ENCODER - The Signal Correlator was designed 
to take these four pressure transducer signals 
and produce one composite continuous pressure 
waveform for the P-V diagram analysis. 

The outputs from the pressure transducer 
are signal-conditioned by Kistler Model 504E 
charge amplifiers. The outputs of the charge 
amplifiers are connected to che four inputs of 
the signal correlator marked PI, P2, P3, and 
P4 (Figure 6). The composite pressure wave- 
form is available at the P-OUT jack. 

An absolute position optical shaft angle 
encoder is used to provide rotor angle infor- 
mation to the signal correlator. The encoder 
is belt -driven from the engine output shaft 
with a 3:2 ratio. Thus, for every three out- 
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put shaft revolutions, the encoder turns two 
revolutions and the rotor, by virtue of inter- 
nal gearing, turns one revoli'tion. The 3 j 2 
encoder drive ratio was chosen so that one 
engine contustion cycle would be represented 
by 720° of encoder rotation as is the case 
for a four cycle piston engine. This allows 
additional NASA-developed instrumentation to 
be used without special modification for this 
engine (12). 

The encoder resolves to 0,35° at its own 
shaft which translates to 0.53° engine out- 
put shaft rotation, and 0.18° rotor rota- 
tion. The active regions and the transition 
(overlap) regions are shown in Figure 7 refer- 
enced to encoder rotation. The 0° encoder 
rotation reference point is chosen such that 
the volume of the combustion chamber being 
investigated is at the minimum volume position 
(non-firing). 

Figure 8 shows (clockwise from upper 
left); (a) the four individual pressure 
transducer signals versus crank angle for one 
cycle where top dead center (TDC) is defined 
as the mini:mum volume condition and bottom 
dead center (BDC) as the maximum volume condi- 
tion so as to be analogous to the piston en- 
gine; (b) the composite pressure signal, 
pieced together and vertically aligned for 
rotor face A for one cycle; (c) the resultant 
pressure-volume diagram displayed on an oscil- 
loscope using the composite pressure signal; 
(d) a 10 times expansion of the P-V diagram to 
better illustrate the pumping loop; and, (e) a 
p-9 blowup of the transducer overlap region. 

Examination of the lOx expansion of 
pumping loop in Figure 8 shows a pressure rise 
approximately midway through the exhaust 
stroke. Further investigation shows that this 
pressure rise began when' the second rotor 
reached the exhaust port. It is believed that 
separate exhaust manifolds would have elimi- 
nated this rise and a 2 % increase in work out 
would have been realized at this speed and 
load. 

The presence of a passage between the 
housing and the pressure transducer diaphragm 
presents an additional complication in inter- 
preting the observed pressure signal. Three 
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such distortions were described by Evers 
(13). They are; phase error due to lap time, 
frictional pressure loss, and a distortion 
caused by response. No attempt was made to 
correct for these errors. Oscillatory pres- 
sures were observed and are shown in Figure 
8E. The question arises as to whether this 
oscillation is caused by a pressure wave in 
the combustion chanter or in the passage of 
the transducer. Using the equation proposed 
by Sears and ^emansky (14) the resonant fre- 
quency, Fr, of this passage length was 
calculated for test point No. 1 as follows; 


L c Characteristic length of passage. 

C s Speed of sound ® temp of the gas. 

The temperature was assumed to be equal to the 
temperature of the diaphragm, 183° F. As- 
suming the gas in the passage is a perfect 
gas, the speed of a wave can be estimated from; 

C = f rgc^f (9) 

® 1800 F ^ C 1243 ft/sec 

. • . Fr = 5736 Hz 

The trace on the photograph (Fig. 8E) 
oscillates at about 5250 Hz which leads one to 
believe that the observed oscillation is in 
the transducer passage. It is also inter- 
esting to note that in this overlap region the 
gas in both transducer passages appears to be 
oscillating at the same frequency. This is to 
be expected as the gas temperature and tube 
length are approximately the same in each of 
the two passages. • 

IMEP INSTRLHENT - The IMEP instrument 
evaluates the j PdV in the crank angle domain 
as follows: 

jL i^PdV - L. / P (iO) 

Vd J Vd 7 l^dey 

c> 

- Crank angle degrees 
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In this the function dV/dB car. he uetpr- 
mineri from the engine georrietry and does not 
need to be computed in real time. This func- 
tion is stored in the Instrument in read-only 
memory (ROM). A plot of dV/dB vs. 8 for the 
engine is shown 3n Figure 9, 

The basic operation of the svstem is il- 
lustrated in the block diagram of Figure 10. 
Engine crank angle data is generated" external 
to the instrument by an absolute shaft angle 
encoder. This crank angle data is input to a 
Memory Address Decoder (MAD). The MAD gener- 
ates the actual memory address f sign of" the 
computation and the computation start and stop 
signals. The 256 valu^is of dV/d8 sta ting at 
bottom dead center (BDC) and ending at top 
dead center (TDC) are preprogrammed and stored 
in Read Only Memory (ROM). Because of the 
symmetry in the volume crank angle relation- 
ship for the rotary engine, the values of 
dV/d0 from TOC to BDC are the negative of 
those from BDC to TDC for equal angles before 
and after TDC. Therefore, only the magnitude 
of dV/d0 is stored in the ROM, This effec- 
tively increases the resolution of a fixed 
memory word size unit by a factor of 2, In 
addition, it allows the MLiltiplying Digital to 
Analog Converter (MDAC) and the Analog to 
Digital Converter (ADC) to be unipolar de- 
vices, thus eliminating any errors due to zero 
effects and zero drift with time or tempera- 
ture. 

The contents of the ROM are applied as the 
digital input to a multiplying or scaling 
digital to analog converter (MDAC). The pres- 
sure signal from' the correlator is applied as 
an analog input to the MDAC. 

The output of the MDAC is an analog signal 
equal to the product P(dV/d8). The signal" is 
then digitized by the ADC. 

The computation of IMEPT is performed 
according to the approximation: 



where K = scaling constant. 

Figure 11 shows a PV diagram for a typical 
computation of p(dV/d8). The computational 
slices contain the following terms: 
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From the intake stroke: V s A1 

Frcm the compression stroke: - P3 V s - (A3 
+ A2 + Al) 

From the power stroke; P4 V s + (AA + A3 
+ A2 + Al ) 

From the exhaust stroke: Po V s - (A2 + 

Al) 

Collecting terms yields the slice sum: A4 - A2 

The terms on the compression and exhaust 
stroke are treated as negative in the computa- 
tion, since (dV/d9) is negative in these re- 
gions. 

The AV si j ces vary in width since the 
computation is performed at constant A9 incre- 
ments. Thus, the slices are more numerous and 
of smaller width in the regions of TDC and 
BDC, This is fortuitous as the region of 
steepest pressure gradients would normally 
occur at these locations. 

In addition to the signal correlator ad- 
ditional instruments were used to measure peak 
pressure and angle of peak pressure. A bar- 
graph generator* was also used to display 100 
consecutive cycles of IMEPT, PMEPp, peak 
pressure and angle of peak pressure for each 
test point. The means and standard deviations 
for these 100 cycle measurements were also 
computed by this instrument for each test 
point. Two typical bargraphs are shown in 
Figure 12A and 12B. Fi^re 12A shows 100 
consecutive cycles obtained in real time using 
the LeRC signal correlator and IMEP meter 
while the engine was running at steady speed 
and load. The upper set of bargraphs on this 
figure represent 100 consecutive cycles of 
IMEPT and the lower represents 100 consecutive 
cycles of PMEP. The second shows the transi- 
tion of IMEPT from firing to motoring. 

RESULTS AND DISCUSSION 

Figures 13, 14, and 1^ are comparisor.s of 
IMEPT calculated from Equation 6 and the same 
quantity itieasured directly. The range of 
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speeds was limited to 30CX) RFM (crankshaft) by 
the motoring capability of our equipment. The 
measurements were taken at the full range of 
loads over which the engine will operate for 
the two speeds tested (see Table 2), 

The average difference between the cyclic 
integral measured and the value of this in- 
tegral back calculated was 6,0%. In addition 
to the fact that we are making measurements on 
only one rotor, there is variability from 
pocket to pocket on the sarre rotor. As we did 
not compute the cyclic integrals of each of 
the three thermodynamic systems simulta- 
neously, this variability is unknown. The 
only manner in which a direct conparison could 
have been made would be to compute the cyclic 
Integral for each of the six systems of inter- 
est simultaneously. For contustion analysis, 
a single rotor engine would he better to use 
than the multirotbr engine as both cost and 
complexity of the necessary instrumentation 
would be reduced. 

Figure 16 shows that agreement between the 
measured and back-calculated values of the 
cyclic integral is poorer at light loads than 
at higher ones. This may not be surprising 
because under light-load conditions small 
differences in rbtor-to-rotor power production 
or system inaccuracies become more pro- 
nounced. In addition, the instrumentation is 
operating at less than 10% of its full scale 
value at" the low-load conditions. 

Figure 17 shows how the average coeffi- 
cient of correlation (Std deviation of IMEPT/ 
Mean IMEPT) varies with power at the two en- 
gine speeds tested. This shows that there is 
more cycle-to-cycle variability at low loads 
than at moderate and high loads. This phenom- 
enon is generally observed in carbureted en- 
gines and is probably due to variations in 
local air fuel ratios. 

Measurements of PMEP were made during 
firing and motoring at constant intake mani- 
fold pressure (IMP), The PMEPp value should 
be exactly equal to PMEP motoring at constant 
IMP and the same speed, if other conditions 
were equal. All other conditions are not 
exactly the same as heat transfer to the in- 
coming air fuel charge will be higher in 
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firinq bhnn Jt\ t-hp tnotorlnq ease* Tills would 
tend to make the averaqe chamber pressure 
during intoKu slightly' higher during firing 
than during motoring and thus reduce PMEP 
during firing. A comparison of PMEP during 
firing and PMEP during motoring is shown in 
Figure 18. 

In general, good agreement between these 
parameters was obtained. The use of a belt- 
driven motoring device conplicated making the 
motoring measurement in that it was difficult 
to maintain speed and constant manifold pres- 
sure, In addition, the engine came equipped 
with a carburetor which included a primary and 
secondary throttle plate. This also made the 
problem of holding manifold vacuum constant 
difficult. 

In Figure IP FMEP is ihown vs. brake power 
for the two nominal speeds tested, FMEP is 
primarily a function of engine speed as would 
be expected. There is more scatter in the 
2000 RPM data than in the 2900 RPM data be- 
cause variation in PMEP has a smaller effect 
in the calculation of FMEP (from the kWEP and 
PMEP measurements) since the frictional term 
becomes more dominant as the speed increases. 


SINMARY AND CXINCLUSIONS 


1. The instrumentation designed to give 
real time IMEP measurement was built and suc- 
cessfully demonstrated. 

2. Cycle -to -cycle variability of the 
cyclic work integral (IMEPT) expressed as std. 
deviation/mean, is greater at low loads than 
at high loads for the rotary engine tested. 

3. An accurate determination of engine 
IMEP will require instrumentation of all the 
systems that produce work. For a two rotor 
engine this will be six simultaneous measure- 
ments of IMEP. 

A. A performance improvement could be 
realized in this engine hy the addition of a 
second exhaust pipe'. 
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PMFP motoring averaged 7% higher than 
PMEP firing for the conditions tested. 

6, The average value of FMep is 31^ 
higher at 2900 RPM than it is at 2000 RPM. 

7. Oscillations in passages do not inter- 
fere with the performance of the correlator in 
a severe manner. 
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Appendix A 

CORRELATOR CIRCUIT DESCRIPTION 

The circuit diagram for the signal correla- 
tor is shown in Figures A1 ami A2. 

Figure A1 shows the analog signal processing 
circuitry. This circuitry contains a repeti-* 
tive circuit block of which the first block is 
composed of Ul, Al, Rl, and Cl. U1 is a 
sample-hold amplifier configured as a unity 
gain inverting surrmer when the sample switch 
is closed. When the sample switch is open, 
the output of Ul is determined by the voltage 
on capacitor Cl. Amplifier Al is configured 
as a unity gain inverting summer. Resistor Rl 
(RIA-RIG) is a laser trimmed thin film resis- 
tor network of seven matched resistors. The 
nominal value of these resistors is lOK + 1^, 
but the seven resistors are matched to o7l%. 
One input to Ul is -PI which is obtained by 
inverting PI with amplifier A5, The other 
input of Ul is P2. Switch SI is closed peri- 
odically during the transition region where 
both PI and P? are reading the combustion 
pressure, ^airing this time the output of Ul 
is the instantaneaus difference between PI and 
P2 (A2. 1) w hen SI is closed, and the average 
value (A2,1) of this difference when SI is 
open. A further discussion of this averaging 
circuit will follow later in this report.' The 
difference (42,1) is applied as one input to 
amplifier Al, The other input of Al is pres- 
sure P2. The output of Al is then: 

-P2^ := P2 - (P2 - PDavG. 

Thus, P2* is P2 offset by the average differ- 
ence existing between P2 and PI during the 
overlap period where both transducers' are 
reading the same pressure, htence, the result 
is to force P2* to equal PI. 

The remaining basic circuit blocks function 
in the same manner. Thus, P3* is forced equal 
to P2* during the overlap region where P3 and 
P2 are both active and PA* is forced equal to 
P3* during their common overlap region. The 
net result is a continuous pressure waveform 
with all transducers referenced back to PI. 
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/Viy offset drift in the individual transducer 
signals is cancelled out during each engine 
cycle. Figure RA shows the four individual 
pressure transducer signals versus encoder 
rotation and Figure 8B shows the composite 
waveform generated by the signal correlator. 

The output of U4 (AX, A) is used as an error 
indicator. Since the transducers are all 
referenced to PI, the final circuit block 
should force PI* to equal PI and A1,A should 
be zero. If this is not the case, then a gain 
or sensitivity error Is present in one or more 
of the transducer waveforms. No attempt has 
been made in the design of the nignal correla- 
tor to correct for gain errors. 

The four signals, -PI, -P2*, -P3*, and -PA*, 
are applied as the inputs to analog multi- 
plexer, U5. This multiplexer selects one of 
the four inputs and passes it to amplifier A7, 
a high input -impedance unity-gain buffer am- 
plifier. Inputs AO and A1 determine which 
input signal is selected. Amplifier A7 is 
used to prevent the multiplexer 'W' resis- 
tances from affecting the output. Amplifier 
A6 is a unity-gain inverting amplifier used to 
restore the polarity of the* input signal. 

Figure A2 shows the circuitry used' to gener- 
ate the control signals for the correlator. 

The major component used for control is an 8 
bit by 256 word read-only memory, U7. The 
address select inputs, AO through A7, are 
driven from the shaft angle encoder. The 
encoder generates digital shaft position sig- 
nals encoded in a Gray code. This Gray code 
is partially converted to a complementary 
natural binary by external circuitry. Addi- 
tional conversion is performed by exclusive 
"or" gates U6A and U6B. The resultant code is 
an 8 bit binary number which specifies the 
encoder shaft angle in increments of 720°/ 

256 or approximately 2.8°. Because the code 
is complementary, encoder position 0° is 
represented by all "one's” and encoder posi- 
tion 720° is represented by all "zeros." 

The outputs of U7 are as follows; 

AO, Al; Used to select which pressure 
signal is passed to the P-OUT 
output . 
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T2»l : Defines the overlap region where 
both transducers PI and P2 are 
active. 

T3,2 : Defines the overlap region where 
both transducers P2 and P3 are 
active. 

TaJ 3 ; Defines the overlap region where 
both transducers P3 and P4 are 
active. 

T1 ,4 ; Defines the overlap region where 

both transducers and PI are 
active. 

The contents of U7 are shown in Table A1 and 
a timing diagram is shown in figure A3. 

Integrated circuits UB and UB are dual one- 
shot multivibrators. The pulse width for all 
four one-shots is set at 10 usee. Bach one- 
shot is triggered from a combination of two 
input signals. One of the inputs is the Tj,i 
signal from the memory and is used to enable 
the one-shot only during the respective over- 
lap region. The other input is BO which is 
generated by the shaft encoder. SO is the 
least significant bit of the encoder position 
and is a square wave signal of period 0.70 
(at 6000 encoder RPM, BB has a period of 2 Dj 
sec). Every time BO goes from logic "0" to 
logic "1" a 10 usee negative going pulse is 
generated by the oneshot if the enabling input 
(Tj,i) is logic "0". These output pulses are 
applied as the sample control inputs to the 
respective sample-hold amplifiers U1 through 
U4. When the control input is logic "0,” the 
sample switch is closed. 

ERROR AVERAGING 


When the switch of the sample-hold amplifier 
is closed, the hold capacitor (Cl, C2, C3, or 
Ca) is charged through the respective series 
resistor. If the switch were closed continu- 
ously, then the voltage across the hold capac- 
itor in response to a step input, E, would 
follow the exponential: 


Equation A1 


Ec = ECo + ( E 

fCo) . 
(l-e“t/T) 
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Where Eq s Voltage across the 

capacitor 

Eqo s Initial voltage across 

the capacitor 

t s Time 

T K Time constant = RC 

However, if the switch were closeii periodi- 
cally only for a brief but fixed period of 
time, t, then the capacitor voltage would 
follow Equation A1 only during that time. At 
all other times the capacitor voltage would 
remain constant at its last value. Therefore, 
the capacitor voltage is a function of the 
number of intervals that tne switch is 
closed. If we let; 

Equation A2 At = T/n 

or 

Equation A3 nAt = T 

n = some fixed number of 

switch closings 

t = period of time the 

switch is closed 

T = time constant 

Then from Equation A3 the tifae constant is 
proportional to the number, n. Since the 
number of switch closures is determined by the 
overlap interval divided by the BO cycle peri- 
od, both being independent of RPM, then the 
filter response or averaging period is RPM 
independent also. 

Having the filter response RPM independent 
allows maxinnum utilization of the overlap 
regions. The hold capacitors are chosen such 
that 5 time constants are contained in each 
overlap interval. 

Thus, Ci is determined from the following 
equctions: 

Equation AA 5n s Tj^/0.7° Schock 

and since; 
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Equation A5 

nAt = T = RCi 

Where 

= number of degrees in a 
specific overlap 
region 

.70 

= period of BO in 
degrees 

n 

= time constant in 
sample periods 

t 

= 10 usee. 

R = 10,000 ohms 

Then substituting Equation A5 into Equation 
AA and solving for Ci yields; 

Equation A6 

Ci = (Ti,j/0.70) 
lAt) (1/5R) 


= Tn (10 X 
10--^)/(0.7O) (5) 
(10 X 103) 


= 28g Tj^i X 10-6 

Ci 

s: Hold capacitor in ufd 

’^j,i 

= Overlap interval in 
encoder degrees 

From the data in 

Table Al, 

T2,1 

= 110.0° Cl = 0.31ufd 


= 17.0 C2 = ,0049ufd 

T4,3 

= 36.0 C3 = .OlOufd 

Ti,A 

=14.0 C4 = .OOAOufd 

The closest standard capacitor values to those 


calculated above were chosen. 
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POWER SUPPLY 


The signal correlator is packaged for con- 
venience’ in a Tektronix TM-500 series blank 
plug-in module kit. The module plugs into a 
Tektronix Power Module (e.g., RTM 506) which 
contains the power transformer, rectifiers, 
filter capacitors, fuses, etc., and also un- 
coirmitted NPN and PNP power transistors for 
use as series pass elements if desired. 

Figure A2 shows the power supply configuration 
used in the signal correlator. Ull is a +15 
volt, 3 terminal regulator, The input is +33 
VDC from the Power Module. C9 is used as a 
filter capacitor. Similarly, U12 is a -15V 
regulator with CIO as a filter capacitor. UlO 
is a +5V regulator which uses one of the ex- 
ternal series pass PNP transistors as the 
series pass element to decrease the current 
flow through the regulator. C12 is an input 
filter capacitor and Cll is an output filter 
capacitor. 01 is a 6.3V zener diode which is 
used to supply -9 VDC to the memory from the 
-15 VDC supply . 

CALIBRATION 

Because of the use of precision matched gain 
resistors, no adjustments are needed for 
gain. The only adjustments are for zero off- 
set. Each sample-hold amplifier has an as- 
sociated potentiometer for offset adjustment 
of both thie sample-hold amplifier and the 
associated summing amplifier. 

The calibration procedure is as follows; 

1. Remove US and U9 from their sockets. 

2. Ground Pin A and Pin 12 of both US and U9 
sockets. This forces all sample-hold switches 
to remain closed. 

3. Apply a variable DC power supply to input 
PI and adjust for -5.000 volts at the output 
of A5 (-P1). 

4. Ground the P2, P3, and PA inputs. 
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5. Adjust for a reading of -5.U’10 volts at 
the output of Al. 

6 . Adjust R9 for a reading of -5.000 volts at 
the output of A2. 

7. Adjust R12 for a reading of -5.000 volts 
at the output of A3. 

8. Adjust R15 for a reading of -5.000 volts 
at the output of AA. 

There is no offset adjustment for A5, A7, or 
A6 because the instrufinents that use the P-OUT 
signal have auto-zero functions incorporated 
in them. 
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Table 1 


Model (MAZDA) 

No. of Rotors 
Displacement inVl 
Rated BMP at RFM 
Compression Ratio 
Ignition 


Engine Description 

1978"12A Non-Bmissions 
? 

70/1.15 
110 P 7000 
9.4:1 

Conventional Breaker Type 



Table 2. Matrix of Test Points 
OTTRW) 


Nominal 

Speed 


2046 

RF#1 

10/7.5 

15/11.2 

30/22.4 

35/26,1 

1, 2, 3, 4 

2046 

RF//2 

10/7.5 

15/11,2 

30/22,4 

35/26.1 

9, 10, 11, 12 

2046 

RF#3 

10/7.5 

15/11.2 

30/22.4 

35/26.1 

17, 18, 19, 20 

2950 

RF//1 

10/7.5 

15/11.2 

30/22.4 

45/33.6 

5, 6, 7, 8 

2950 

RF//2 

10/7,5 

15/11.2 

30/22.4 

45/33.6 

13, 34, 15, 16 

2950 

RF#5 

10/7.5 

15/11.2 

30/22.4 

45/33.6 

21, 22, 23, 24 


*RF = Rotor Face 



TABLE A1 - ROM DATA 


Encoder Degrees 
(To Nearest Whole Degree) 

00 

I 

3 

6 

I 

17 

20 

I 

115 

118 

I 

225 

228 

i 

329 

332 

I 

346 

349 

I 

515 

518 

I 

551 

554 

I 

717 


HEX Address 

l-EX 

FF 

F3 

1 

I 

FE 

F3 

FD 

1 

73 

i 

1 

F9 

1 

73 

F8 

1 

FO 

1 

1 

D6 

1 

FO 

D5 

1 

EO 

1 

1 

AF 

1 

EO 

AE 

1 

FI 

1 

8A 

1 

FI 

89 

( 

01 

1 

84 

1 

Dl 

83 

1 

F2 

j 

48 

I 

F2 

47 

1 

B2 

1 

1 

3B 

B2 

3A 

1 

F3 

1 

1 

00 

1 

F3 



Figure AL, - Correlator schematic lanalogl 














Rgure A2. - Correlator schematic (digital plus Pwr supl. 
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Fi('ureA3. « Timing diagram. 



Figure lA. - Transducer locations. 
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Figure IB. - Transducer active regions. 
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Figure 5. - Transducer inslallallon in 
rotor housing. 
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Figure 6. * Photograph of signal correlator. 
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Figure 7. - Transducer active regions. 
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Figure 8. - Transducer and signal correlator outputs. Engine condition ?OOC RPM, 38 ft*lbf, Hi • .88. 












Figure 9. - Rate of change of volume with crank angle. 
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Figure 14. >• Back calculation of IMEPT and measured IMEPT 
versus power, 
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Figure 15» - Back calculation of IMEPT and measured IMEPT 
versus power. 
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Figure 16. - Average error in IMEPT measurements, 







